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Abstract 

A Voigt profile function emerges in several physical investigations 
(e.g. atmospheric radiative transfer, astrophysical spectroscopy, plasma 
waves and acoustics) and it turns out to be the convolution of the Gaus- 
sian and the Lorentzian densities. Its relation with a number of special 
functions has been widely derived in literature starting from its Fourier 
type integral representation. The main aim of the present paper is to 
introduce the Mellin-Barnes integral representation as a useful tool to 
obtain new analytical results. Here, starting from the Mellin-Barnes 
integral representation, the Voigt function is expressed in terms of the 
Fox H-function which includes representations in terms of the Mei- 
jer G-function and previously well-known representations with other 
special functions. 

PACS numbers: 32.70.-n, 02.30.Gp, 02.30.Uu 



1 Introduction 

During the recent past, representations of the Voigt profile function in terms 
of special functions have been discussed in this Journal [1, 2, 3, 4, 5, 6]. 
With the present paper we aim to continue those researches. In particular, 
we derive the representation in terms of Fox H-function, which includes 
previously well-known results. 
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The Voigt profile function emerges in several physical investigations as 
atmospheric radiative transfer, astrophysical spectroscopy, plasma waves 
and acoustics and molecular spectroscopy in general. Mathematically, it 
turns out to be the convolution of the Gaussian and the Lorentzian densities. 
Here we are studying the ordinary Voigt function and not its mathematical 
generalizations, for example, see the papers [7, 8, 9]. The computation of 
the Voigt profile is an old issue in literature and many efforts are directed 
to evaluate this function with different techniques. In fact, an analytical 
explicit representation in terms of elementary functions does not exist and 
it can be considered a special function itself. Moreover it is strictly related 
to the plasma dispersion function [10] and to a number of special func- 
tions as, for example, the confluent hypergeometric function, the complex 
complementary error function, the Dawson function, the parabolic cylinder 
function and the Whittaker function, see e.g. [1, 2, 3, 11, 12, 13, 14]. All 
previous representations are derived starting from the integral formula due 
to Reiche in 1913 [15] that is actually a Fourier type integral. 

The Voigt profile function remains nowadays a mathematically and com- 
putationally interesting problem because computing profiles with high ac- 
curacy is still an expensive task. The actual interest on this topic is proven 
by several recent papers (2007) on mathematical [16, 17, 18, 7, 19, 20] and 
numerical aspects [21, 22, 23]. A huge collection exists of published works 
on this topic, but instead to report it, we give the significative datum that 
searching in Google Scholar the strings "voigt profile function" and "voigt 
function" the number of files found is ~ 24, 600 and ~ 70, 100, respectively. 

The Mellin-Barnes integrals are a family of integrals in the complex plane 
whose integrand is given by the ratio of products of Gamma functions. De- 
spite of the name, the Mellin-Barnes integrals were initially studied in 1888 
by the Italian mathematician S. Pincherle [24, 25] in a couple of papers on 
the duality principle between linear differential equations and linear differ- 
ence equations with rational coefficients. The Mellin-Barnes integrals are 
strongly related with the Mellin transform, in particular with the inverse 
transformation. As shown by O.I. Marichev [26], the problem to evaluate 
integrals can be successfully faced with a powerful method mainly based on 
their reduction to functions whose Mellin transform is the ratio of product 
of Gamma functions and then, after the inversion, the problem consists in 
the evaluation of Mellin-Barnes integrals. Moreover, they are also the es- 
sential tools for treating higher transcendental functions as Fox H-function 
and Meijer G-function and a useful representation to compute asymptotic 
behaviour of functions [27]. 

The main object of the present paper is to introduce the Mellin-Barnes 
integral representation as a useful tool to obtain new analytical results that 
in the future can lead to efficient numerical algorithms for the Voigt function. 
A successful application of such approach has been shown in [28], where the 
parametric evolution equation of the Voigt function (and its probabilistic 
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generalization) is derived and the scaling laws, with respect to the parame- 
ter, in the asymptotic regimes are computed using the Mellin-Barnes integral 
representation. After all, this work can be seen also as an interesting exercise 
per se of application of the Mellin-Barnes integral method. 

The rest of the paper is organized as follows. In section 2 the basic def- 
inition of the Voigt profile function is given and some classical and recent 
representations are reviewed. In section 3 the Mellin-Barnes integral repre- 
sentation of the Voigt function is derived and, in section 4 starting from this 
result, first the Voigt function is expressed in terms of the Fox H-function 
and later, in cascade, the representation with the Meijer G-function and 
other special functions are obtained. Finally in section 5 the summary and 
conclusions are given. 



2 The Voigt profile function 

2.1 Basic definition 

The Gaussian G(x) and the Lorentzian L(x) profiles are defined as 



Gix) = —= — exp 



L(x) = ji^hr, (l) 



where ujq and ujl are the corresponding widths. The variable x is a wave- 
number and then its physical dimension is a length raised up to —1. The 
Voigt profile V{x) is obtained by the convolution of G{x) and L(x) 

V(x) = L(x- 0G(0 da = -^jf / 2 dt . (2) 

The comparison between the Voigtian, Gaussian and Lorentzian profile is 
shown in figure 1 for different values of the width ratio ujl/ujg- 
Let /(k) be the Fourier transform of f(x) so that 

/+oo -1 /*+oc 

^e +i ™f(x)dx, f{x) = —J^ e-^f(K)dK, (3) 

then 

V(k) = G(k)L(k) = e -^ 2 /4 e -^|K| ? ( 4) 

and 



1 f + OO 

V(x) = — / e -«^ e -^« 2 /4~^l«l dK 
1 r+oo 

= - e -^ L K-U> G K /4 cos ( KX ) d ^ (5) 

TT Jo 
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a=0.01 



a=0.1 




Figure 1: Comparison between the Voigtian, Gaussian and Lorentzian pro- 
files with a = lu l /u; g = 0.01, 0.1, 1, 2. 
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Formula (5) is the integral representation due to Reiche [15]. 

2.2 Some classical and recent representations 

Let x 1 be the dimensionless variable x' = x/ujq, the Voigt function can be 
re-arranged in the form 

V(x) = -^K(xf, y) , K(x, y) = V -\ - 6 2 <£ , (6) 

where y = ujl/ujg and from (3) it follows that 

1 f +°° 

K(x, y) = - r e-yt-t l A cos(xO . (7) 

V 71 " Jo 

The Voigt function does not possess an explicit representation in terms of 
elementary functions and several alternatives to (2) have been given in lit- 
erature, mainly with the intention to obtain a more efficient numerical com- 
putation. 

Combining x and y in the complex variable z = x — iy, the function 
K(x,y) (6) is 

• r+oc -£ 2 

K(x,y) = Re[W(z)}, W{z) = - -d£. (8) 

7T J_ oa Z~i 

where is strongly related to the plasma dispersion function [10] and 

some classical representations can be found, see e.g. [11, 12, 13]. In fact, 
form the relation of W(z) with the complex complementary error function 
Erfc(-iz) and the Dawson function F(z) = e~ z J Q Z e^ d£ then it follows 
that 

K(x,y) = Re[W{z)], W(z)=e~ z2 Erfc(-iz), y>0, (9) 

K(x, y) = Re[W{z)\ , W(z) = e~ z ' 2 + ^=F{z) . (10) 

More recent representations are, for example, those derived in 2001 by 
Di Rocco et al [29] 



*<*•»> -B- 1 )" {ferny 1 * ( 

n=0 K \ ' J \ 



2n + l 1 2 



2 '2' 
2a . / .3 



1F1 (n + l,-;y 2 \}x 2n , (11) 
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where \F\{a,(3\ z) is the confluent hypergeometric function, and in 2007 by 
Zaghloul [20], which "completed" a previous formula given by Roston & 
Obaid [30], 



K(x, y) = [1 - erf (y)\ e^W) cos ( 2 xy) + f e^ x2+ ^ sm[2y(x - £)] d£ . 

V 71 " Jo 

(12) 

Recently (2007), He & Zhang [17] claimed to have derived an exact 
calculation of the Voigt profile that is proportional to the product of an 
exponential and a cosine function. However this representation assumes 
negative values in contrast with the non-negativity of the Voigt function. 
For this reason, this result is not correct. 

Further representations are given in terms of special functions, see for 
example the one involving the confluent hypergeometric function \F\ [2, 3] 

K{x,y) = e^ 2 -* 2 )cos(2xy) - {(y + ix) 1^(1; 3/2; (y + ix) 2 ) + 

(y-ix) 1 F 1 (l;3/2;(y-ix) 2 )} , (13) 

and others involving the Whittaker function Wk t m, the Erfc— function and 
the parabolic cylinder function [14, formulae (17,13,16)] 

K(X, y) = ^L={(y- ix)' 1 / 2 e^— )V2 W _ 1/4 ,_ 1/4 ((y - ixf) + 

(y + ix)' 1 / 2 e (f+») 2 / 2 W_ 1/4i _ 1/4 ((y + ix) 2 )} , (14) 
K(x, V) = \ [e (y ~ ix)2 Erfc(y - ix) + e (y+ix)2 Erfc(y + ix)} , (15) 



K(x, y) = 6 {e~ lxy D^[V2(y - ix)] + e +lxy D_i[v^(y + »*)]} 

\/'2ir I. J 



2vr 

(16) 

3 The Mellin-Barnes integral representation 

Let us consider again dimensional variables and the Gauss, Lorentz and 
Voigt functions defined as in (1-2). From (4) we have 

V(x) = — / e -^ e -^« 2 /4-^|K| dK = _ | /+ ( x ) + ; ( 17 ) 

2vr J^oo 27T 
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where 

r+oo 

I±(X)= e -^L±ix)K-^l± dK _ 

JO 

Following [27], the Mellin-Barnes integral representations of I±(x) can be 
obtained from the definition of the Gamma function 

r+oo 

r(z)= e 2 "V«de, 

J 

and the Mellin-Barnes integral representation of e~ z (z / 0) 

e~ z = — [ r(s)z- s ds = Y tl^z n , 

where C denotes a loop in the complex s plane which encircles the poles 
of T(s) (in the positive sense) with endpoints at infinity in Re(s) < and 
with no restrictions on argz [27]. The functions I±(x) have the following 
Mellin-Barnes integral representations 



I±(X) = ^ +0Oe "^ K2/4 {^- / r(s)[{LJ L ±ix)K]- S dsjdK 
= r ( S ){^ +0 ° e "^ K2/4K " Sds } (LO L ±ix)- S dK 



i i 

WG 27TZ J C 



(uj l ± ix) 



(18) 



Hence the Mellin-Barnes integral representation of the Voigt function is 
given by 



V(x) 



1 



2ttloq 2iri J £ 



— (uj l +ix) 



ds + 



1 

2tH 



T(s)T 



1 s 

2 ~ 2 



uj g 



■{ujl - ix) 



ds 



(19) 



However we note that for a complex number z = \z\e ie , 9 = arct&n{I m{z) / Re{z)), 
the following rule holds 



z n + z n = \z\ n e ine + \z\ n e- ine 
= \z\ n {e ind + e~ in9 ) 

= 2\z\ n cos(n8) = 2|z| n cos(narctan(Im(2;)/i?e(z))) , 
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and formula (19) becomes 



V{x) 



1 1 



ttujg 2m J £ 



— r(,)r 



cos 



s arctan 



-s/2 



(is . 



G 



(20) 

Consider again (18), changing s — > — s and taking the corresponding 
integration path £ as a loop in the complex plane that encircles the poles of 
r(— s), an other Merlin-Barnes integral representation of the Voigt function 
equivalent to (19) is 



V(x) = 



2m J c y 1 \2 ^2 
2m J c K ; V2 2 



— (ujl + ix) 



LOG 



-(uj l - ix) 



ds + 



ds 



(21) 



and in more compact form 



V(x) = — -L / T(-s)T (\ + 

V ' TTUJG 2lTlJ C V ' \2 



COS 



s arctan 



x 



,uj t + ar 



s/2 



ds . 



G 



(22) 



4 The Fox H- and Meijer G-function representa- 
tions 

4.1 The Fox H-function representation 

The Voigt function can be represented in terms of well known special func- 
tion, see §2.2, but its representation in terms of Fox H-function is still not 
known. The expression in terms of Fox H-function is important because 
it is the most modern representation method and, actually, it is the most 
compact form to represent higher transcendental functions. The definition 
of the Fox H-function is given in Appendix. 

For ujg and ujl fixed, 2(u>l + ix)/u>G / and 2(u>l — ix)/u>G / 0, from 
(21) and (19) we have, respectively, 



V(x) = 



2ttujg\ 11 


2 / 

— (lo l + IX) 
WG 


(1/2,1/2) " 
(0,1) 


+ 






— (uj l - IX) 


(1/2,1/2) ■ 
(0,1) . 


}■ 


(23) 
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V(x) 



2ttujg 



H 



W 11 



_2{u L + ix) 
_2(lu l - ix) 



(1,1) 
(1/2,1/2) 

(1,1) 
(1/2,1/2) 



+ 



(24) 



As a consequence of the fact that this is the most comprehensive represen- 
tation, in cascade, the others with less general functions can be obtained. 

4.2 The Meijer G-function representation 

The first H-function in (23) can be rewritten in terms of the Meijer G- 
function. Setting Z = 2(ujl + ix)/log, we have 



H 



(1/2,1/2) 
(0,1) 



l/ £ r(- 2s) r(I + s )^ 

V^27riJ c V ; V2 
Z 2 



s ! I I - + .s 



—= Lx 

V7T 



12 



1/2 
0,1/2 



Z< 



2\ s 



(25) 



where the change of variable s — > 2s and the duplication rule for the Gamma 
function, T{2z) = T{z)T{l/2 + z)2 2z - 1 tt~ 1 / 2 , are applied. The second H- 
function in (23) follows from the first with Z = 2{lol — ix)/uJG- Finally, the 
Voigt function in terms of the Meijer G-function is given by 



V(x) = 



2vr 3 / 2 



n 2\ 
^12 



1/2 
0,1/2 



Z 2 



1/2 
0,1/2 



(26) 



The Meijer G-function can be reduced to other special functions and, for 
example, representations given in (14, 15, 16) are straightforwardly recov- 
ered. 

In fact, the Meijer G-function and the Whittaker function W^^m are 
related by the formula [31, p. 435] 



= T(b - a + l)r(c - a + l)x ib+c - l V 2 e x ' 2 W Km {k) , 

where k = a — (b + c + l)/2 and m = (b — c)/2. Thus the G-functions in 
(26) are expressed in terms of the Whittaker function as 



^21 
^12 





a 


z 


b, c 
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n 2\ 
^12 



Z 2 1/2 
0,1/2 



Z 



'1/2 



e Z2/8 ^_ 1/4 ,_ 1/4 (Z 2 /4). 



Hence, the Voigt function in terms of the Whittaker function is 

-1/2 



V(x) = 



1 Z 



e Z %-i/4,-i/4(2 2 /4)+ 



-1/2 



Moreover, from the identity 



G 21 



12 



\z 2 


1/2 


_ 4 


0,1/2 



2/ e^/ 8 ^ 1/4 ,_ 1/4 (Z 2 /4)^. (27) 



7re z2 / 4 Erfc(Z/2), 



we have 



= [e^/'ErfciZ^ + e^^ErfciZ/l)} , (28) 



and from 



we have 



^12 



rz 2 


1/2 


_ 4 


0,1/2 



27re z / 8 D-^Z/y/2), 



V(x) = -^J— [e z2 / 8 D^(Z/V2) + e 22 / 8 D^Z/V^)} . (29) 
Setting log = 1 in (27,28,29) formulae (14,15,16) are recovered, respectively. 

5 Summary and conclusions 

In the present paper the Mellin-Barnes integral representation of the Voigt 
profile function is derived. We think that this integral representation is 
a useful tool to have new analytical and numerical results in the subject. 
Starting from this, the Voigt function has been expressed in terms of the Fox 
H-function, which is the most comprehensive representation. In cascade, the 
expression in terms of the Meijer G-function is obtained and the previous 
well-known representations with the Whittaker, the Erfc and the Parabolic 
cylinder functions are recovered. 
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Appendix. The Fox H- and the Meijer G- functions 

The Fox H-function is defined as [32, 33] 



TTran 
n pq 



(a 1 ,A 1 ), ■ ■ ■ , (a p ,A p ) 
(b 1 ,B 1 ),---,(b g ,B q ) 



2m 



h(s)z s ds , 



where 



h(s) 



ut=i - Bjs) n u r(i - aj + A jS ) 



U q j= m + i r(i - bj + b jS ) u p J=n+ i r K' - A > s ) ' 

an empty product is interpreted as unity, {m , n , p , q} are nonnegative 
integers so that < m < q, < n < p, {Aj , Bj} are positive numbers and 
{aj , bj} complex numbers. The H-function is an analytic function of z and 
makes sense if [32, 33] 

i) ^ , n > 

a) o < |z| < p- 1 , // = o 

where 



p 9 
.A, 



j=i j=\ j=i j=i 

The Meijer G-function corresponds to the special case Aj 



Bk 



(j = 1, ...,p;k= l,...,q) 



Tiran 

pq 



, (dp, 1) 



(6i,l),---,(6„l) 



— ^pq 



ai , ■ ■ ■ , a p 
h,---,b q 
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